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We present a precision measurement of the top-quark mass using the full sample of Tevatron 
^/s = 1.96 TeV proton-antiproton collisions collected by the CDF II detector, corresponding to 
an integrated luminosity of 8.7 fb~^. Using a sample of ti candidate events decaying into the 
lepton+jets channel, we obtain distributions of the top-quark masses and the invariant mass of 
two jets from the W boson decays from data. We then compare these distributions to templates 
derived from signal and background samples to extract the top-quark mass and the energy scale 
of the calorimeter jets with in situ calibration. The likelihood fit of the templates from signal and 
background events to the data yields the single most-precise measurement of the top-quark mass, 
Mtop = 172.85 ± 0.71 (stat) ± 0.85 (syst) GeV/c^ 

PACS numbers: 14.65. Ha, 13.85.Ni, 13.85.Qk, 12.15.Ff 



The top quark (t) is by far the heaviest known el- 
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ementary particle [ij. It contributes significantly to 
electroweak radiative corrections relating the top-quark 
mass (Mtop) and the W boson mass to the mass of the 
Higgs boson 0. Precision measurements of Aftop pro- 
vide therefore important constraints on the Higgs boson 
mass. Over the last decade, the CDF and DO experiments 
have been improving the precision of the Mtop measure- 
ment [3] , joined recently by the experiments at the Large 
Hadron Colhder [|. 

This Letter reports the single most-precise measure- 
ment of the top-quark mass to date using the template 
method In this measurement, we reconstruct top- 

quark masses in each event and compare the distribution 
of data with template distributions derived from model 
calculations to estimate Mtop. We also use the template 
distributions of hadronically-decaying W bosons to con- 
strain the jet energy scale, an important uncertainty in 
the Mtop measurement. This is an update of a previous 
measurement that used a subset of the present data and 
determined Mtop =172. 2± 1.5 GeV/c^ H- This measure- 
ment has an accuracy approximately 30% better than an 
earlier result with similar template technique and 12% 
better with respect to the previous best measurement . 
In this measurement, we use not only larger samples but 
also improve jet energy calibration using an artificial neu- 
ral network [10] to achieve better jet energy resolution, 
and increase signal acceptance allowing less pure signal 
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samples into the analysis. The measurement is performed 
with ^/s = 1.96 TeV Tevatron proton-antiproton collision 
data collected by the CDF II detector [ll| and corre- 
sponding to an integrated luminosity of 8.7 fb~^. This is 
the full Run II CDF data set with requirements of good 
performances in the all detector elements. 

Assuming unitarity of the quark- mixing matrix [l| , the 
top quark decays almost exclusively into a W boson and 
a b quark. In the ti events, the case where one W decays 
leptonically into an electron (e) or a muon {^) plus a 
neutrino (i^), including the cascade decay oi W ^ tv 
and T — > fiv or r — > ev, and the other W decays 
hadronically into a pair of jets (tt — > lubbqq) defines 
the lepton-|-jets decay channel. Lepton-l-jets events are 
selected by requiring one isolated jl2l | electron (muon) 
with Et > 20 GeV {pt > 20 GeV/c) and pseudorapidity 
|?7| < 1. 1 [ 13l| . We also require large missing transverse 
energy > 20 GeV, and at least four jets. Jets 

are reconstructed with a cone algorithm [ist with radius 
R = v^(A?7)2 + (A(/))2 = 0.4. In addition to the stan- 
dard jet energy corrections Jl(ij, we train a neural net- 
work including additional information to the calorimeter 
one, such as jet momentum from the tracker as described 
in Ref. (T^]. We have performed the training separately 
for b quarks and light flavor quarks from W boson de- 
cay. The additional information on the jets improves 
the resolution of the reconstructed top-quark masses and 
hadronically decaying W boson mass. This allows ap- 
proximately 13% better statistical precision compared 
with standard correction. Jets originating from b quarks 
are identified (tagged) using a secondary- vertex tagging 
algorithm (T^]. We divide the sample of candidate lep- 
ton-l-jets events into sub-samples based on the number of 
identified b jets, zero 6-tagged jet (0-tag), one 6-tagged 
jet (1-tag), and two or more 6-tagged jets (2-tag). In the 

0- tag events, we require exactly four tight jets (transverse 
energy Et > 20 GeV and \r]\ < 2.0). In the 1-tag and 
2-tag events, three tight jets and at least one loose jet 
(i?T > 12 GeV and \r]\ < 2.4) are required. We divide 
the 1-tag and 2-tag samples into sub-samples based on 
the number of tight jets and call the "tight" sub-sample 
the one requiring exactly four tight jets and "loose" the 
one consisting of the remaining events passing selection. 
The measurement uses five sub-samples (0-tag, 1-tagL, 

1- tagT, 2-tagL, and 2-tagT, where L and T represent 
loose and tight selection, respectively). The sub-samples 
of 0-tag and 1-tagL are newly added in this measure- 
ment. The introduction of these higher-background sub- 
samples offers a gain of approximately 12% in the sta- 
tistical precision. We apply an additional requirement 
on the scalar sum of transverse energies in the event, 
Ht = i;^°P*™ + ^T + Ejots-^T* El' to further reject 
backgrounds, where e}^^^°'^ is electron or muon trans- 
verse energy or momentum, respectively, and Ei^^ is the 
transverse energy of the jet. 

The primary background sources are W-l-jets and QCD 
multijet production. We also consider small contribu- 



tions from Z-l-jets, diboson, and single-top-quark pro- 
duction. To estimate the contribution of each process, 
we use a combination of data- and Monte Carlo (MC)- 
based techniques described in Ref. fl^. For the Z-fjets, 
diboson, single top, and events we normalize simulated 
events using their theoretical cross sections [20l-|2a. We 
use the data-driven techniques described in Ref. (23| to 
estimate QCD multijet background. The shape of the 
VF-|-jets background is obtained from simulation, while 
the number of H^-|-jets events is determined from the 
data by subtracting the other backgrounds and tt event 
contributions. 

For each event, three observables are used, two recon- 
structed top-quark masses (m™'^° and mj;"™'^'') and the 
invariant mass of the two jets from the hadronically de- 
caying W boson (mjj). We have a complete reconstruc- 
tion of the ti kinematics in the lepton-|-jets channel Q 
with constraints from the precisely known W boson mass 
and requiring the t and t masses to be the same. Assum- 
ing that the leading four jets in the detector originate 
from the ti decay products, there are twelve, six, and 
two assignments of jets to quarks for 0-tag, 1-tag, and 
2-tag events, respectively. A minimization is performed 
for each assignment using a comparison to the ti hy- 
pothesis with ml°'^° taken from the assignment that yields 
the lowest (Xmin)- To increase the statistical power 
of the measurement, we employ an additional observable 
^icco(2) assignment that yields the second low- 

est x^. To reject poorly reconstructed events, we require 
Xmin < 3 and Xmin < 9 for and tagged (both 1- 

tag and 2-tag) events, respectively. The dijet mass TOjj is 
calculated as the invariant mass of two non 6-tagged jets 
that provides the closest value to the known W boson 
mass of 80.39 GeV/c^ 24]. We apply boundary con- 
ditions on mi'=™ and ml'"'"'^^^ (100 < mf'^™, m^''™^^^ < 
350 GeV/c^) and also mjj (60 < mjj < 110 GeV/c^, 
50 < TOjj < 120 GeV/c2, and 50 < mjj < 125 GeV/c^ 
for 0-tag, 1-tag, and 2-tag events, respectively). The es- 
timated number of background events and the observed 
numbers of events after event selection, x^, and boundary 
requirements are listed in Table HI 

We estimate the probability density functions (p.d.f.s) 
of signal and background using a kernel density estima- 
tion method [25| . A three-dimensional kernel density 
estimation Q accounts for the correlation between the 
three observables. The dijet mass tojj of the two jets as- 
signed to the W in the lepton-|-jets channel is used for 
in situ calibration of the jet energy scale (JES) [1, Q. 
The p.d.f.s for the observables are estimated at 76 dis- 
crete values of Mfp from 130 GeV/c^ to 220 GeV/c^ 
and at 29 discrete values of Ajes from —3.0 CTc to 3.0 ctc, 
where CTc is the CDF JES fractional uncertainty based 
on a combination of instrumental calibration and analy- 
sis of data control samples [11] . The parameter Ajes is 
defined as the difference between the true jet energy and 
our measured jet energy after calibration by a factor of 
1 -I- Ajes • 0, @| ■ We interpolate the MC distributions 
to find p.d.f.s for arbitrary values of Mtop and Ajes us- 
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TABLE I: Expected and observed numbers of signal and background events assuming tt production cross-section a^t = 7.45 pb 
and Mtop = 172.5 GeV/c^ 





0-tag 


1-tagL 


1-tagT 


2-tagL 


2-tagT 


W+jets 


703 ± 199 


170 ± 60 


102 ± 37 


11.6 ± 4.9 


8.4 ± 3.5 


Z+jets 


52.3 ± 4.4 


8.9 ± 1.1 


5.9 ± 0.7 


0.8 ± 0.1 


0.5 ± 0.1 


Single top 


4. a ± v.o 


iU.i3 ± U.9 


U.O ± U.D 


on 1 no 
Z.Z ± U.O 


1.1 ± U.z 


Diboson 


60.3 ± 5.6 


111 ± 1.4 


8.5 ± 1.1 


1.0 ± 0.2 


0.8 ± 0.1 


Multijcts 


143 ± 114 


34.5 ± 12.6 


20.7 ± 16.6 


4.4 ± 2.5 


2.5 ± 2.4 


Background 


963 ± 229 


235 ± 61 


144 ± 41 


19.9 ± 5.5 


13.8 ± 4.2 


tt signal 


645 ± 86 


695 ± 87 


867 ± 108 


192 ± 30 


304 ± 47 


Expected 


1608 ± 245 


930 ± 106 


1011 ± 115 


212 ± 30 


318 ± 47 


Observed 


1627 


882 


997 


208 


275 



ing the local polynomial smoothing method 26]. We fit 
the signal and background p.d.f.s to the distributions of 
the observables in the data using an unbinned maximum 
likelihood fit [13, . Independent likelihoods are used 
for each sub-sample, 0-tag, 1-tagL, 1-tagT, 2-tagL, and 2- 
tagT, and the total likelihood is obtained by multiplying 
them together ^-Q. References provide detailed 

information about this technique. 

We test the mass determination using 1 500 statistical 
trials for a set of 11 different Mtop values ranging from 
160 GeV/c^ to 185 GeV/c^. In each experiment, we draw 
the number of signal and background events each from 
a Poisson distribution centered at the expected number 
of signal and total background shown in Table IH respec- 
tively. The distributions of the average deviation from 
the input top mass and the width of the deviation nor- 
malized to the estimated uncertainty for simulated exper- 
iments are corrected to be unity and zero, respectively. 
The correction is M^°" = 1.03 x M™<='^" - 4.88 GeV/c^, 
where M^'^'^ is the maximum likelihood estimate and 
M^°" is the corrected value of the measurement. We in- 
crease the measured uncertainty by 2.9% to correct the 
width of the pull. We also test the mass fit results using 
different values of Ajes between —1.0 dc to 1.0 ac with 
three different Mtop points, 168, 173, and 178 GeV/c^. 
With a correction discussed above, the residuals of Mtop 
from different Ajes values are consistent with zero. 

We examine various sources of systematic uncertain- 
ties by comparing the analysis results of statistical trials 
in which we vary relevant parameters within their uncer- 
tainties. The dominant sources are the residual JES [l^ 
and signal modeling. We vary the JES parameters within 
their uncertainties in both signal and background MC 
generated events and interpret the deviations of the re- 
sults as additional uncertainties. The uncertainty aris- 
ing from the choice of MC generator (signal modeling) 
is estimated by comparing the results of pseudoexper- 
iments generated with pythia and HERWIG [30i] . 
We examine the effects of higher-order corrections us- 
ing MC@NLO [sTl , a full next-to-leading-order simulation, 
which is added in the signal modeling. The system- 



TABLE 11: Estimated systematic uncertainties (units in 
GeV/c^). 



Source 


Systematic uncertainty 


Residual jet energy scale 


0.52 


Signal modeling 


0.57 


6 jet energy scale 


0.18 


b tagging efficiency 


0.03 


Initial and final state radiation 


0.06 


Parton distribution functions 


0.08 


Gluon fusion fraction 


0.03 


Lepton energy scale 


0.03 


Background shape 


0.20 


Multiple hadron interaction 


0.07 


Color reconnection 


0.21 


MC statistics 


0.05 


Total systematic uncertainty 


0.85 



atic uncertainty to the energy corrections of b jets (6- 
JES) arising from our modeling of b fragmentation, b 
hadron branching fractions, and calorimeter response, 
captures the additional uncertainty not taken into ac- 
count in the residual JES. The uncertainty on the b- 
tagging efficiency can propagate into a bias of the Mtop 
measurement which is taken as a systematic uncertainty. 
The uncertainty due to the limited knowledge of initial- 
state radiation is constrained by studies of radiation in 
Drell-Yan events. We vary both initial and final state ra- 
diation within these constraints by extrapolating in the 
Pt of Drell-Yan events to the tt mass region Q. We 
estimate the systematic uncertainty due to parton dis- 
tribution functions by varying the independent eigenvec- 
tors of the CTEQ6M [s^] parton distribution functions, 
varying the QCD scale (Aqcd), and comparing our nom- 
inal CTEQ5L [H with MRST72 0] parton distribution 
functions. We vary the gluon fusion fraction from 5% to 
20%, corresponding to the one standard deviation upper 
bound on the gluon fusion fraction 35]. We estimate sys- 
tematic uncertainties due to the lepton energy and mo- 
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FIG. 1: Distributions of the three variables ml'"^°, mjj, and ml'"^°^^\ used to measure Mtop for 0-tag and tagged events. The 
data are overlaid with the predictions from the kernel density estimation probability distributions using Mtop = 173 GeV/c^ 
and the full background model. 



mcntum scales by propagating shifts in electron energy 
and muon momentum scales within their uncertainties. 
The background shape systematic uncertainty accounts 
for the variation of the background composition. We es- 
timate a multiple-hadron-interaction systematic uncer- 
tainty to account for the difference in the average number 
of interactions between simulation and data. The color 
reconnection systematic uncertainty [36| is evaluated us- 
ing samples simulated with and without color reconnec- 
tion effects with different pythia tunes [37*] . All system- 
atic uncertainties are summarized in Table |TT1 The total 
systematic uncertainty adding individual components in 
quadrature is 0.85 GeV/c^. The details of systematic 
uncertainty evaluations are in Ref. ^5,, jjj . 

We perform the likelihood fit to the data and apply the 
corrections obtained using the simulated experiments, 
and measure 

Mtop = 172.85 ±0.71 (stat) ± 0.85 (syst) GeV/c^. 

Figure [1] shows the distributions of the observables used 
for the Mtop measurement in the lepton-t-jets channel, 
overlaid with density estimates using tt signal events with 
-Mtop = 173 GeV/c^ (close to the measured Mtop) and the 
full background model. 

In conclusion, we have performed a measurement of 
the top-quark mass using the template method in the 
lepton-|-jets using the full CDF Run II data set corre- 
sponding to 8.7 fb~^ pp collisions. The result, Mtop = 



172.85 ± 1.11 GeV/c^, is the best single measurement of 
this important physics parameter. It is consistent with 
the most recent Tevatron average of Mtop = 173.18 ± 
0.94 GeV/c2 [| and will significantly contribute to the 
Tevatron's and world's average value for the top quark 
mass. 
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